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Abstract 



This paper presents measurements of the polarization of W bosons in top quark decays, derived 
from tt events with missing transverse momentum, one charged lepton and at least four jets, or two 
charged leptons and at least two jets. Data from pp collisions at a centre-of-mass energy of 7 TeV 
were collected with the ATLAS experiment at the LHC and correspond to an integrated luminosity of 
1 .04 fb _1 . The measured fractions of longitudinally, left- and right-handed polarization are F = 0.67± 
0.07, F L = 0.32 ± 0.04 and F R = 0.01 ± 0.05, in agreement with the Standard Model predictions. As 
the polarization of the W bosons in top quark decays is sensitive to the Wtb vertex Lorentz structure 
and couplings, the measurements were used to set limits on anomalous contributions to the Wtb 
couplings. 
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1 Introduction 

Discovered in 1995 by the CDF and D0 experiments [1, 2], the top quark is the heaviest 
known fundamental particle, with a mass of 173.2±0.9 GeV [3]. Measurements of top quark 
properties play an important role in testing the Standard Model and its possible extensions. 
One particular test is the study of the Wtb vertex Lorentz structure and couplings, which 
can be probed by measuring the polarization of W bosons produced in top quark decays. 

At the LHC, top quarks are produced mainly in pairs via the strong interaction and 
are predicted to decay via the electroweak interaction into a W boson and a bottom quark 
with a nearly 100% branching fraction. Events with ti pairs can thus be classified according 
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to the decay of the two produced W bosons. Each boson can decay either into a quark- 
antiquark pair or into a charged lepton and a neutrino. The single-lepton and dilepton 
topologies, both considered in the analyses presented in this paper, have one and two 
isolated charged leptons in the final state. Only electrons and muons, including those from 
r decays, are considered here. 

The Wtb vertex is defined by the electroweak interaction and has a (V — A) structure 
where V and A are the vector and axial-vector contributions to the vertex. Since the W 
bosons are produced as real particles in top quark decays, their polarization can be longitu- 
dinal, left-handed or right-handed. The fractions of events with a particular polarization, 
Fq, Fl and Fr, are referred to as helicity fractions. They are predicted in next-to- next- 
to-leading-order (NNLO) QCD calculations to be F = 0.687 ± 0.005, F L = 0.311 ± 0.005, 
Fr = 0.0017 ± 0.0001 [4]. These fractions can be extracted from measurements of the 
angular distribution of the decay products of the top quark. The angle 6* is defined as 
the angle between the momentum direction of the charged lepton from the decay of the 
W boson and the reversed momentum direction of the 6-quark from the decay of the top 
quark, both boosted into the W boson rest frame [5]. The angular distribution is: 

= 7 i 1 ~ cos2 r ) F ° + I 0- ~ cosr ) 2 F L + | (1 + cos A*) 2 F R . (1.1) 
a d cos 0* 4 '8 8 

All previous measurements of the helicity fractions, performed by the CDF and D0 Col- 
laborations [6-8] at the Tevatron, are in agreement with Standard Model predictions. 

Information about the polarization of the W bosons can also be obtained through 
complementary observables, such as the angular asymmetries, A + and defined as: 

N(cobO* >z)- N{cos6* < z) 
± ~ iV(cosfl* > z) + N{cos6* < z) ' 

with z = ±(1 — 2 2 / 3 ) for A±, allowing the dependence on F^ and Fr to cancel, respectively. 
The asymmetries can be related to the helicity fractions by a simple system of equations [9, 
10]. In the Standard Model, the NNLO values for these asymmetries are A + = 0.537±0.004 
and A_ = -0.841 ± 0.006 [4]. 

In the presence of anomalous Wtb couplings the helicity fractions and angular asymme- 
tries depart from their Standard Model values [5, 10]. In effective field theories, dimension- 
six operators can be introduced which modify the Wtb vertex [11-13]. Coefficients con- 
trolling the strength of these operators can be constrained by measurements of the helicity 
fractions or the angular asymmetries. 

This paper describes measurements of the W boson polarization in top quark decays 
and the constraints on the Wtb vertex structure based on a data set recorded with the 
ATLAS detector between March and June 2011 and corresponding to an integrated lu- 
minosity of 1.04 fb _1 . The helicity fractions were measured using two different methods. 
The first compares the observed cos#* distribution with templates for different W boson 
helicity states obtained from simulation. The second method extracts angular asymme- 
tries from an unfolded cos#* spectrum corrected for background contributions. Limits on 
anomalous couplings, generated by the aforementioned dimension-six operators, were set 
using the combined result from the two measurements. 
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This paper is structured as follows. Data samples and simulated samples modelling 
signal and background processes are described in Section 2, followed by a summary of the 
event selection in Section 3. The analysis strategies can be found in Section 4. Sources of 
systematic uncertainty are discussed in Section 5, the results are summarized in Section 6 
and Section 7 concludes the paper. 

2 Data and simulation samples 

The ATLAS detector [14] at the LHC covers nearly the entire solid angle around the colli- 
sion point. It consists of an inner tracking detector surrounded by a thin superconducting 
solenoid, electromagnetic and hadronic calorimeters, and an external muon spectrometer 
incorporating three large superconducting toroid magnet assemblies. A three-level trigger 
system, designed to reduce the event rate from 40 MHz to about 200 Hz, is used to select 
events of interest. 

2.1 Data sample 

Data from pp collisions at a centre-of-mass energy of 7 TeV were collected using single- 
lepton triggers with transverse momentum thresholds of 20 GeV for electrons and 18 GeV 
for muons. The size of the data sample corresponds to an integrated luminosity of 1.04 ± 
0.04 fb" 1 [15, 16]. 

2.2 Signal and background modelling 

Signal and most background processes were modelled by Monte Carlo (MC) simulation. 
The signal process was simulated using two different generators. The leading order (LO) 
Monte Carlo generator Protos [10, 17] was used with the CTEQ6L1 sets of parton dis- 
tribution functions (PDFs) [18] to generate three samples assuming Fo = 1, = 1 and 
Fr = 1. This was achieved by choosing appropriate values for anomalous Wtb couplings, 
as described in Refs. [9, 10]. The output was interfaced to Pythia [19] to simulate par- 
ton showers and hadronization. In addition, the next-to-leading-order (NLO) generator 
MC@NLO [20-22] was used with the CTEQ6.6 PDF set [23] for studies of systematic 
differences in the top quark production and decay modelling. 

The W and Z boson production in association with multiple jets was simulated using 
the Alpgen generator [24] and the CTEQ6L1 PDF set. Tree-level matrix elements with 
up to five final-state partons were included. The standard Alpgen matching scheme 
was used to remove overlaps between the n and n + 1 parton samples. Heavy flavour 
samples containing Wbb, Wcc, Wc and Zbb events were simulated separately. The Z+jets 
samples were generated with dileptons in the invariant mass range 10 < mn < 2000 GeV. 
Diboson processes were simulated using Herwig [25]. Single top quark production was 
simulated using MC@NLO, invoking the 'diagram removal scheme' [26] to remove overlaps 
between the single top quark and tt final states. An additional sample of VF+jets events 
was generated with Sherpa [27] and used to study systematic uncertainties. Details are 
described in Section 5. Apart from the Protos tt samples, all events were hadronized with 
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Herwig using Jimmy [28] for the underlying event model. Different underlying event tunes 
were used, depending on the hadronisation program used (i.e. Pythia or Herwig) [29]. 

Additional M C samples were used for the evaluation of systematic uncertainties. These 
were generated with the AcerMC [30] and Powheg [31, 32] generators interfaced to 
Pythia for hadronization. They are described in Section 5. A top quark mass of 172.5 GeV 
was assumed for all signal samples if not stated otherwise. Additional MC@NLO samples 
were generated assuming different top quark masses. 

All simulated events contain multiple pp-interactions and contain pileup contributions 
corresponding to a bunch spacing of 50 ns. These simulated events were re-weighted such 
that the average number of interactions per proton-proton bunch crossing was the same in 
data and MC simulation. The average number of interactions per event was roughly six. All 
samples were processed by the detector and trigger simulation after event generation [33, 
34], and subjected to the same reconstruction algorithms as the data. 

The cross-section of simulated tt samples was normalized to 164.6 pb, the value ob- 
tained from approximate NNLO calculations [35-37], and consistent with recent measure- 
ments [38]. While the Z+jets normalization was obtained from the NNLO QCD cross- 
section calculations in the single lepton channels, a data-driven estimate for the normaliza- 
tion was used for the dilepton channels [38] . The diboson with jets production was rescaled 
to match NLO calculations of the inclusive cross-sections. 

For the single-lepton analysis, the multijet production background, where an electron 
or muon originates from hadron decay or instrumental background, was estimated from 
data as described in Refs. [39, 40]. The VF+jets background was obtained from simulation, 
except that its normalization was derived from data, from a study of the asymmetry in the 
production of W + and W~ bosons [40]. For the dilepton analysis, background contributions 
from F^+jets, single-lepton ti and single top quark production were estimated using the 
method described in Ref. [38]. In the single-lepton and dilepton analyses this source of 
background is labelled as "misidentified leptons" . 

3 Event selection 

3.1 Physics object definition 

The reconstruction and identification of electrons, muons, jets and the magnitude of the 
missing transverse momentum (i?™ lss ) used in the analyses presented here followed the 
criteria employed for the measurement of the ti production cross-section [39]. More details, 
including a description of the trigger and the vertex requirements, can be found in Ref. [40]. 

The ti events with at least one isolated charged lepton (electron or muon) in the final 
state were considered as signal. Events with r leptons decaying into muons or electrons in 
the final state were also considered to be part of the signal. 

All the considered events were required to fulfil general event quality criteria and to 
have a well-defined primary vertex with at least five associated tracks. Reconstructed 
electrons were required to have £r > 25 GeV and \rj\ < 2.47 excluding 1.37 < \r/\ < 1.52, 
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while muons were required to have pr > 20 GeV and |r/| < 2.5 . Jets were reconstructed 
with the anti-fci algorithm [41] with a radius parameter equal to 0.4. They were required 
to have pi > 25 GeV and |r/| < 2.5. For electrons, the energy not associated to the electron 
cluster but contained in a cone of AR = a/(A0) 2 + (A??) 2 = 0.2 around it was required 
not to exceed 3.5 GeV. For muons, the sum of track transverse momenta and the total 
energy deposited in a cone of AR = 0.3 around the muon were both required to be less 
than 4 GeV. Muons reconstructed within a AR = 0.4 cone of a jet satisfying > 20 GeV 
were removed to reduce the contamination caused by muons from hadron decays within 
jets. Subsequently, jets within AR = 0.2 of an electron candidate were removed to avoid 
double counting, which can occur because electron clusters are also reconstructed as jets. 

A combination of two algorithms was used for 6-jet identification (6-tagging). They are 
based on reconstructed secondary vertices and the impact parameter significances of tracks 
within jets. The chosen working point resulted in a 6-tagging efficiency of 70% for jets 
originating from 6-quarks in a sample of simulated tt events and a light-quark jet rejection 
factor of about 100 [42]. 

3.2 Single-lepton channels 

The single-lepton tt signal is characterized by a high-px isolated charged lepton and missing 
transverse momentum from the neutrino, from the leptonically decaying W boson, two 
light-quark jets from the hadronically decaying W boson and two 6-quark jets. The two 
channels with either an electron or a muon in the final state are referred to as single-electron 
and single-muon channels. The following event selection requirements were applied: 

• the appropriate single-electron or single-muon trigger had fired; 

• events were required to contain exactly one isolated electron or muon; 

• in the single-electron channel, E!f lss > 35 GeV and m^iW) > 25 GeV were required 2 
while in the single-muon channel the criteria were E™ 1SS > 20 GeV and E™ ss + 
m T (W) > 60 GeV; 

• events were required to have at least four jets, with at least one of them being tagged 
as a 6-jet. 

The numbers of events expected after all selection requirements are shown in Table 1 
for the single-electron and single-muon channels, together with the observed numbers of 
events in the data sample. The single-electron channel contains significantly fewer events 
than the single-muon channel due to the more stringent requirements on the lepton E^ 

1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the 2-axis along the beam pipe. The :r-axis points from the IP to the centre 
of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, cj>) are used in the transverse 
plane, <f> being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 9 as r\ = — lntan(#/2). 

2 Here mxfW) is the W boson transverse mass, reconstructed as ^/2pi r pY r [1 — cos((// — 0")] where the 
measured missing transverse momentum provides the neutrino information. 
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Table 1. Event yields in the single-electron and single- muon channels after the event selection. 
The table shows the expected number of events including their estimated total uncertainty as well 
as the number of events observed in the data sample. 



Process 


Single electron 


Single muon 


tt 


4400 ± 1100 


6500 ± 1400 


T^+jets 


900 ± 700 


1400 ± 1000 


Z+jets 


120 ± 90 


140 ± 90 


Diboson 


14 ± 12 


22 ± 12 


Single top 


260 ± 90 


360 ± 110 


Misidentified leptons 


220 ± 220 


500 ± 500 


Total predicted 


5900 ± 1300 


9000 ± 1800 


Data 


5830 


9121 



(j)t), E!f lss and m^(W). The observed numbers of events are in agreement with those 
expected, with purities of 75% and 72% for the single electron and single muon channels. 
The systematic uncertainty discussed in Section 5 dominates the error on the event yields, 
with jet energy scale, signal and background modelling, and 6-tagging being the dominant 
sources of uncertainty. Example distributions in data are compared with the Standard 
Model predictions in Figure 1. The data are described well by the expected distributions. 

3.3 Dilepton channels 

The dilepton ti signal is characterized by two high-px isolated charged leptons, missing 
transverse momentum corresponding to the undetected neutrinos from the two leptonically 
decaying W bosons, and two 6-quark jets. The three channels are referred to as ee, \x[i and 
efj, channels. The following event selection requirements were applied: 

• the single-electron trigger had fired for the ee channel, the single-muon trigger had 
fired for the fifi channel and either of these triggers had fired for the e/x, channel; 

• events were required to contain exactly two oppositely charged and isolated leptons 

(ee, fifj, or e/i); 

• at least two jets, with at least one of them being 6-tagged, were required; 

• in order to avoid the low-mass Drell-Yan background region, events were required to 
have mu > 15 GeV; 

• in the ee and /i/i channels, the missing transverse momentum had to satisfy E™ ss > 
40 GeV, and the invariant mass of the two leptons had to differ by at least 10 GeV 
from the Z boson mass, i.e. \mu — mz\ > 10 GeV, with mz = 91 GeV, in order to 
suppress backgrounds from Z+jets events and events containing misidentified leptons; 
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Figure 1. The top row shows (left) the Er spectrum of the electron in the single-electron channel 
and (right) the pt spectrum of the muon in the single-muon channel. The bottom row shows (left) 
the missing transverse momentum in the single-electron channel and (right) the pt spectrum of the 
leading jet in the single-muon channel. The error band shows the statistical uncertainty from the 
Monte Carlo simulation. The last bin in each distribution includes the overflow. 

• in the e/i, channel, no E™ ss or Z boson mass veto cuts were applied; however, the 
event H^, defined as the scalar sum of the transverse energies of the two charged 
leptons and all selected jets, must satisfy > 130 GeV to suppress backgrounds 
from Z+jets production; 

• events were required to pass the ti reconstruction criteria, as described in Section 4.1. 

The numbers of events expected after all selection requirements are shown in Table 2 
for the three dilepton channels, together with the observed numbers of events in the data 
sample. The observed numbers of events is in agreement with those expected, with a 
purity of 94% for all three channels. The systematic uncertainty discussed in Section 5 
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Table 2. Event yields in the dilepton channels after the event selection and event reconstruction. 
The table shows the expected number of events including their estimated total uncertainty as well 
as the number of events observed in the data sample. 



Process 


ee channel 


fifi channel 


efx channel 


tt 


159 ± 28 


320 ± 40 


750 ± 100 


Z+]ets (ee,/i//) 


0.7 ± 2.0 


2 ± 4 




Z+jets (tt) 


0.0 ± 3.8 


1.4 ± 1.3 


4.7 ± 3.1 


Diboson 


0.2 ± 0.3 


0.6 ± 0.3 


1.3 ± 0.5 


Single top 


5.0 ± 1.3 


10.8 ± 2.4 


24 ± 5 


Misidentified leptons 


4 ± 2 


7 ± 4 


19 ± 9 


Total predicted 


170 ± 30 


340 ± 40 


800 ± 100 


Data 


191 


354 


836 



dominates the error on the event yields, with jet energy scale, jet energy resolution, signal 
and background modelling, and 6-tagging being the dominant sources of uncertainty. Ex- 
ample distributions in data are compared with the Standard Model predictions in Figure 2. 
The data are described well by the expected distributions. 

4 Measurement strategies 

The two methods to extract the helicity fractions are described below. The analysis was 
carried out for the single-electron, single-muon, ee, efi and \i\x channels separately. The 
reconstruction of tt events in the single-lepton channels differs for the two methods: the 
template method relies more strongly on the structure of the cos 9* distribution than the 
calculation of the angular asymmetries. It thus relies on an improved estimate of the 
particle energies. A common reconstruction method was used in the dilepton channels. 
The results of the analyses and their combination are described in Section 6. 

4.1 Measurement using the template method 

Templates for different signal and background processes were fitted to the observed cos 9* 
distributions based on events reconstructed with a kinematic fit in the single-lepton chan- 
nels [40, 43]. The fit was based on a likelihood which took into account the Breit-Wigner 
forms of the lineshapes of the top quark and the W boson as well as the energy resolution 
of the measurements of the jets and the charged leptons. The missing transverse momen- 
tum was identified with the x- and y-components of the neutrino momentum. The top 
quark mass was fixed to 172.5 GeV, and the mass of the W boson was fixed to 80.4 GeV. 
The likelihood was maximized with respect to the energies of the final state quarks for 
each association of jets to quarks. The permutations were weighted according to the b- 
tagging information and the weights were derived from the efficiency and mis-tag rate of 
the 6-tagger. The permutation with the largest value of the likelihood was used in the 
following steps of the analysis. The overall efficiency for the reconstruction of the correct 
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Figure 2. The top row shows (left) the Et spectrum of the leading electron in the ee channel and 
(right) the pt spectrum of the leading muon in the channel. The bottom row shows (left) the 
missing transverse momentum in the efi channel and (right) the pt spectrum of the leading jet in 
the efi channel. The error band shows the statistical uncertainty from the Monte Carlo simulation. 
The last bin in each distribution includes the overflow. 

event topology was found to be 74% in simulated tt events. Only those events where four 
jets and a lepton are matched to partonic particles, corresponding to roughly 30% of the 
events, were considered for the efficiency computation. 

For the dilepton channels, the tt system was reconstructed by solving a set of six 
independent equations: 

n Ul 4- n U2 — R miss - n^l 4- rF* — Rmiss. /_ i _ \2 _ 2 . 
Px +Px — i Py +Py — ' KPli+Pui) — m Wi 

{Pe 2 + Pvif = m w\ (PW! + Ph) 2 = mf; {pw 2 + p h ) 2 = m\. (4.1) 

I?™ 1SS and E™ 1SS represent the x- and y-components of the missing transverse momentum, 
•pn x and pg 2 (pj 1 and pj 2 , p Ul and p U2 ) correspond to the four- momenta of the two charged 
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leptons (jets, neutrinos). Here, mw and m t are the W boson and top quark masses, 
respectively. For events with only one 6-tagged jet, this jet and the non-tagged one with 
the highest were taken as 6-jets from the top quark decays. If at least two 6-tagged 
jets are present, the two 6-tagged jets with the largest pt in the event were used. The 
pairing of the jets with the charged leptons was based on the minimization of the sum of 
the invariant masses rri£ 1 j 1 and m,£ 2 j 2 . Simulations show that this criterion gives the correct 
pairing in 68% of the events. Up to four solutions can be found to the Equations 4.1. The 
solution with the minimum product of neutrino transverse momenta was chosen and is 
motivated by the low px-spectrum of the neutrinos. In case no solution was found, the top 
quark mass parameter was varied in the 157.5-187.5 GeV range, to try to find one or more 
solutions. In case a range of top masses provided a solution, that with the top quark mass 
parameter closest to 172.5 GeV was taken. Again, multiple ambiguities were resolved by 
taking the solution with the minimum product of neutrino transverse momenta. If still no 
solution was found, the second pairing of jets and charged leptons was used. If no solution 
was found after this procedure the event was discarded. This happened in about 25% of 
the events in data and in simulated ti events. 

In the single-lepton channels, the following six processes were used as templates: three 
different helicity state signal processes, the background from events with misidentified 
leptons, the VT+jets contribution and the sum of all other sources of background. The 
three signal templates were derived from simulation and include single-lepton and dilepton 
events which pass the event selection. These distributions are shown in Figure 3. The 
W+jets template was obtained from simulation. The template for events with misidentified 
leptons was obtained from data. Other processes contributing to the observed spectrum 
include single top quark, diboson and Z+jets production and were summed in a single 
template obtained from simulation. 

In the dilepton channels, templates for the three different helicity state signal processes 
were used as well as a single template representing the different background contributions, 
namely single top quark production, processes with misidentified leptons, production of 
Z bosons with additional jets and diboson production. All templates were obtained from 
simulation, except the misidentified leptons template, which was obtained from data. 

A binned likelihood fit was used to estimate the expected number of events contributing 
to the distribution, assuming independent Poisson-distributed fluctuations in each bin. The 
number of events expected in the ith bin, Aj, is given as: 

AW 

A,= £ A.^+]TA/, (4.2) 
h=-l,0,+l j=i 

where the first term describes the expected number of signal events at particle level with 
defined helicity h. These are reduced by reconstruction efficiencies and acceptances summa- 
rized by the factors e h estimated from simulation. The second term describes the expected 
number of events from background processes and iVbkg is the number of background tem- 
plates. The contributions of the various background templates were constrained according 
to the central values and uncertainties given in Table 1 and Table 2. The uncertainties 
were assumed to be uncorrelated. 
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Figure 3. Distributions of cos 6** for each of the three simulated signal templates. The templates 
for the combined (left) single-lepton and (right) dilepton channels are shown. 

4.2 Measurement using the angular asymmetries 

The angular asymmetries were measured using the cos 9* distribution obtained by recon- 
structing the events with a x 2 minimization technique in the single-lepton channels [44]. 
The x 2 was defined according to: 



2 _ {™ivja ~ m t) 2 {m jbjcjd -m t ) 2 (m £u - m w f {m jcjd - m w ) 2 

X o ' 2 ' 2 ' 2 ' {^^J 

where m t = 172.5 GeV, raw = 80.4 GeV, at = 14 GeV and aw = 10 GeV are the expected 
top quark and W boson mass resolutions, £ represents the selected electron or muon, m^ u 
is the invariant mass of the electron (muon) ctnd the neutrino, ctnd jab,c d 

corresponds to all 

possible combinations of four jets among all selected jets in the event (with m£ U j a , rrij b j c j d 
and rnivja being the corresponding invariant masses). The neutrino was reconstructed using 
the missing transverse energy, with the longitudinal component of the neutrino momentum 
(p v z ) allowed to vary. The solution corresponding to the minimum x 2 value is chosen. No 
6-tagging information was used in resolving the ambiguities in assigning jets. 

The method described in the previous subsection was used to reconstruct dilepton tt 
events. 

For the measurement of the angular asymmetries, A + and A_, the cos 6* distribution 
was divided into four non-uniform bins, which were used to count the number of events 
above and below z = ±(1 — 2 2 / 3 ), as defined in Equation 1.2. A background subtraction 
in the observed cos 9* distribution was performed. Subsequently, the following steps were 
applied iteratively: the number of reconstructed events above and below cos 9* = z were 
counted in data for each asymmetry and correction factors were evaluated by comparing 
the Standard Model expectation with the reconstructed number of simulated ti events. 
These factors allowed corrections to be made for event selection and reconstruction effects. 
The obtained angular asymmetries were then converted into W boson helicity fractions 
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and these values were used to re-derive the cos 6* distribution and evaluate new correction 
factors. The procedure was repeated until the method converged, i.e. until the differences 
between the output observables and the input hypothesis at the previous iteration step 
were below 0.5%. Closure tests were performed using Monte Carlo samples with different 
helicity fractions and no bias was observed. 

5 Systematic uncertainties 

5.1 Sources of systematic uncertainty 

Several sources of systematic uncertainty were taken into account in the analyses presented 
here. These were categorized into the modelling of the signal and background processes, 
and the detector modelling. The impact of the systematic uncertainties on the final results 
is summarized in Table 3. 

5.1.1 Signal and background modelling 

The signal process was modelled with different Monte Carlo generators. Sources of sys- 
tematic uncertainty considered here were the choice of generator and parton shower model, 
the choice of parton distribution functions, the assumed top quark mass and the choice 
of parameters which control the amount of initial and final state radiation. Predictions 
from the MC@NLO and Powheg generators were compared. The parton showering was 
tested by comparing two Powheg samples interfaced to Herwig and Pythia, respec- 
tively. The amount of initial and final state radiation was varied by modifying parameters 
in AcerMC interfaced to Pythia. The parameters were varied in a range comparable 
with those used in the Perugia Soft/Hard tune variations [45]. The impact of the choice 
of parton distribution functions was studied using the reweighting procedure described in 
Ref. [44]. MC@NLO samples were generated assuming different top quark masses and 
their predictions were compared. The observed differences in the results were scaled to 
variations of 0.9 GeV in the top quark mass according to the uncertainty on its Tevatron 
average value [3]. The impact of different models of colour reconnection was studied by 
comparing samples simulated with AcerMC using the Perugia 2010 tune with and without 
colour reconnection [45] as well as the tune A-Pro and ACR-Pro [46, 47]. 

Background processes were either modelled by simulation or were estimated in aux- 
iliary measurements. The number of events with misidentified leptons was estimated in 
data for each channel [40] and the uncertainty on the normalization was estimated to 
be 50% before and 100% after the 6-tagging requirement. The normalization of W+jets 
processes was estimated from supplementary measurements using the asymmetric produc- 
tion of positively and negatively charged W bosons. The uncertainty was estimated using 
Berends-Giele-scaling [48] which yielded 48% for events with four jets and increased with 
the jet multiplicity by 24% per additional jet [39]. Systematic uncertainties on the shapes 
of the T^+jets distributions were assigned based on samples with different simulation pa- 
rameters such as the minimum transverse momentum of the parton and the functional 
form of the factorization scale in Alpgen. Scaling factors correcting the fraction of heavy 
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flavour contributions in simulated VF+jets samples were estimated in auxiliary measure- 
ments described in Ref. [40]. The uncertainties were 76% for Wbb+jets and W^cc+jets 
contributions, and 35% for Wc+jets contributions. The uncertainty on the normalization 
of Z+jets events was estimated using Berends-Giele-scaling. The uncertainties in the nor- 
malization were 48% for events with four jets and increased with the jet multiplicity by 
24% per additional jet. A systematic uncertainty in the shape was accounted for by com- 
paring simulated samples generated with Alpgen and Sherpa. The uncertainty on the 
normalization of the small background contributions from single top quark and diboson 
production was estimated to be about 10% (depending on the channel) and 5%, respec- 
tively. The former estimate was based on the difference between the predictions from 
MC@NLO and MCFM [49], whereas the latter comes from scale and PDF uncertainties 
evaluated with MCFM. 

For the template method, the bin content in each template was varied according to 
a Poisson distribution to estimate the impact of the finite Monte Carlo sample size used. 
For the measurement of the angular asymmetries, the tt reconstruction parameters were 
varied. In the single-lepton channels, the top quark and W boson mass resolutions used in 
the x 2 definition were changed by 25%. For the dilepton channels, the assumed top quark 
mass window used in the kinematic equations was reduced by 50%, taking the difference 
in the observables with respect to the nominal procedure as the uncertainty estimate. The 
uncertainty due to the finite Monte Carlo sample size was found to be negligible since fewer 
bins were used in the unfolding than for the templates. These sources of uncertainty are 
labelled "method-specific uncertainties" in Table 3. 

5.1.2 Detector modelling 

The mis- modelling of lepton trigger, reconstruction and selection efficiencies in simulation 
was corrected for by scale factors derived from auxiliary measurements of the processes Z — > 
and Z — > ee [50, 51]. The uncertainties were evaluated by changing the event selection of 
the supplementary measurement and by testing the stability of the results against changing 
LHC and ATLAS run conditions. The same processes were used to measure the lepton 
momentum scale and resolution. Scale factors and their uncertainties were derived to match 
the simulation to observed distributions. Details are given in Ref. [39]. 

The jet energy scale was derived using information from test-beam data, LHC collision 
data and simulation. Its uncertainty varies between 2.5% and 7% in the central r\ region, 
depending on jet pp and rj [52] . This includes uncertainties in the flavour composition of the 
samples and mis-measurements from close-by jets. An additional pp-dependent uncertainty 
of up to 2.5% was assigned to jets matched to 6-quarks (using Monte Carlo generator- 
level information) due to differences between light-quark and gluon jets as opposed to jets 
containing fo-hadrons. Additional uncertainties of up to 5% (8%) in the central (forward) 
region were added due to pileup. The energy resolution for jets in Monte Carlo simulation 
was adjusted to that observed in data. Uncertainties on the energy resolution of 4-45%, 
decreasing with jet pt, were assigned. The reconstruction efficiency of jets was found to 
be in good agreement with the predictions from simulation and uncertainties of 1-2% were 
assigned. 
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Table 3. Sources of systematic uncertainty and their impact on the measured W boson helicity 
fractions for the combined single-lepton and dilepton channels. The systematic uncertainties were 
symmetrized by using the larger uncertainty. 



Source 


Uncertainties 






F 


F L 


Fr 


Signal and background modelling 


Generator choice 


0.012 


0.009 


0.004 


ISR/FSR 


0.015 


0.008 


0.007 


PDF 


0.011 


0.006 


0.006 


Top quark mass 


0.016 


0.009 


0.008 


Misidentified leptons 


0.020 


0.013 


0.007 


W+jets 


0.016 


0.008 


0.008 


Other backgrounds 


0.006 


0.003 


0.003 


Method-specific uncertainties 


0.031 


0.016 


0.035 


Detector modelling 


Lepton reconstruction 


0.013 


0.006 


0.007 


Jet energy scale 


0.026 


0.014 


0.012 


Jet reconstruction 


0.012 


0.005 


0.007 


6-tagging 


0.007 


0.003 


0.004 


Calorimeter readout 


0.009 


0.005 


0.004 


Luminosity and pileup 


0.009 


0.004 


0.005 


Total systematic uncertainty 


0.06 


0.03 


0.04 



The uncertainties on the momenta of electrons, muons and jets were propagated into 
the missing transverse momentum. A 10% uncertainty from pileup was added in addition. 

The 6-tagging efficiencies and mis-tag rates have been measured in data [42]. Jet p^- 
dependent scale factors, applied to simulation to match the data, have uncertainties which 
range from 9% to 16% and 12% to 45%, respectively. 

The uncertainty on the measured luminosity was estimated to be 3.7% [15]. 

Due to a hardware failure, a small, rectangular region of the ATLAS electromagnetic 
calorimeter could not be read out in a subset of the data (0.87 fb" 1 ). Data and Monte 
Carlo events in which a jet or an electron were close to the affected calorimeter region 
were rejected. The systematic uncertainty labelled "calorimeter readout" in Table 3 was 
evaluated by varying the criteria to reject these events. 

6 Results 

The two methods discussed in Section 4 were applied to the data set described in Section 2. 
Figure 4 shows the observed distribution of cos 8* in the single-lepton and dilepton channels 
together with the sum of the templates scaled to the best fit parameters obtained from the 



-14- 




Figure 4. Distributions of the reconstructed cos 9* used in the template method for data (markers), 
fitted background (dotted line), the Standard Model prediction (dashed line) and the best fit value 
(solid line) for the (left) singlc-lcpton and (right) dilcpton channels. The total uncertainties on the 
helicity fractions for the best fit values are represented by the grey band. For the dilepton channels, 
each event contributed with two entries, corresponding to the two leptonic decays of the W bosons. 



template method as well as the Standard Model expectations. Figure 5 shows the distribu- 
tion of cos 9* in the single-lepton and dilepton channels after the background subtraction 
and the correction for detector and reconstruction effects as used by the asymmetry method 
as well as the Standard Model expectations. 

The combination of the individual measurements of the W boson helicity fractions (Fq 
and Fl) and asymmetries (^4+ and A^) in the single-lepton and dilepton channels was 
done using the best linear unbiased estimator (BLUE) method [53, 54]. The impact of the 
systematic uncertainties described in Section 5 was studied and the results are summarized 
in Table 4, which also gives the combination of all channels for each method. The results 
are compatible with each other and with the final combination. In addition, the template fit 
was repeated with Fr fixed to zero 3 . With the precision of the current measurements, this 
differs negligibly from the Standard Model value, and also follows the approach suggested 
in Ref [13]. The results can also be found in Table 4. 

The results for the angular asymmetries from the single-lepton channels are A + = 
0.52 ± 0.02 (stat.) ± 0.03 (syst.) and A_ = -0.84 ± 0.01 (stat.) ± 0.02 (syst.) whereas the 
results for the dilepton channels are A + = 0.56 ± 0.02 (stat.) ± 0.04 (syst.) and A_ = 
-0.84 ± 0.02 (stat.) ± 0.04 (syst.). 

Most of the measurements and the combined result are limited by systematic uncer- 
tainties. The largest sources of uncertainty are the signal and background modelling, as 
well as the jet energy scale and jet reconstruction. The template fits are more sensitive 

3 In the evaluation of the angular asymmetries no assumption is made for the helicity fractions, so it is 
not possible to fix Fr to zero. 
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Figure 5. Unfolded distributions of cos 8* for the (left) single- lepton and (right) dilepton channels. 
The error bars on unfolded data (markers) include both the statistical and systematic contributions. 
For comparison, the Standard Model NNLO QCD prediction (dashed line) and its uncertainty [4] 
are also shown. 



Table 4. Summary of the W boson helicity fractions measured using the two different techniques 
described and the combination. The quoted uncertainties are the statistical (first) and the system- 
atic (second) uncertainties. 



Channel 


F 


F L 


Fr 


W boson helicity fractions from the template fit 


Single leptons 


0.57 ±0.06 ±0.09 


0.37 ±0.03 ±0.04 


0.07 ±0.03 ±0.06 


Dileptons 


0.92 ±0.10 ±0.10 


0.17 ±0.06 ±0.07 


-0.09 ± 0.05 ± 0.06 


Combination 


0.66 ± 0.06 ± 0.07 


0.33 ±0.03 ±0.03 


0.01 ± 0.03 ± 0.06 


i*R fixed 


0.66 ± 0.03 ± 0.04 


0.34 ±0.03 ±0.04 


(fixed) 


W boson helicity fractions from the angular asymmetries 


Single leptons 


0.66 ±0.03 ±0.08 


0.33 ± 0.02 ± 0.05 


0.01 ± 0.01 ± 0.04 


Dileptons 


0.74 ±0.06 ±0.10 


0.27 ±0.03 ±0.05 


-0.01 ±0.03 ±0.05 


Combination 


0.67 ±0.04 ±0.07 


0.32 ± 0.02 ± 0.04 


0.01 ± 0.02 ± 0.04 


Overall combination 


0.67 ±0.03 ±0.06 


0.32 ± 0.02 ± 0.03 


0.01 ± 0.01 ± 0.04 



to shape-related uncertainties, such as ISR/FSR and jet reconstruction, while the angu- 
lar asymmetries are more sensitive to background normalization uncertainties, due to the 
background subtraction which needs to be performed. 

6.1 Combination 

The results presented in the previous sections were combined using the BLUE method. 
Both the statistical correlations between analyses, and the correlations of systematic uncer- 
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tainties, were taken into account. The statistical correlations were estimated to be between 
40% and 49% for the single-lepton channels and between 83% and 89% for the dilepton 
channels using pseudo-experiments obtained from simulated samples. The background- 
related systematic uncertainties were assumed to be fully correlated within single-lepton 
channels and within dilepton channels, but uncorrelated between single-lepton and dilep- 
ton measurements. The only exception is the uncertainty due to misidentified lepton back- 
ground which depends on the lepton flavour and thus belongs to the group of lepton-related 
uncertainties, which were assumed to be fully correlated between the channels with same 
flavour leptons. The method systematic uncertainties were assumed to be uncorrelated 
between channels and the remaining sources of systematic uncertainty were assumed to be 
fully correlated between channels. Various tests were performed in which the correlations 
among the sources of systematic uncertainty were varied. It was found that the values 
assumed for the correlations were conservative. The systematic uncertainties on the com- 
bined values are summarized in Table 3. The four measurements of the helicity fractions 
and the combined values are shown in Table 4 and Figure 6. 

The individual measurements agree reasonably well within their total uncertainties. 
The x 2 /dof for the global combination of the template fit and asymmetries measurements 
was 0.8 with a % 2 -probability of 75%, where dof is the number of degrees of freedom. The 
largest difference between two measurements is that between the single-lepton and dilepton 
channels obtained with the template method. Since the measurements were performed in 
five independent channels (single electron, single muon, ee, efj, and fifi), the combination 
was performed based on the five individual measurements taking into account all correla- 
tions. The x 2 /dof calculated using the BLUE method for this combination was 1.3 with 
a x 2 -probability of 23%. 

The combined W boson helicity fractions are: 

F = 0.67 ± 0.03 (stat.) ± 0.06 (syst.) , 
F L = 0.32 ± 0.02 (stat.) ± 0.03 (syst.) , 
F R = 0.01 ± 0.01 (stat.) ± 0.04 (syst.) . 

The correlation coefficient between Fq and i*L was estimated to be —0.96. For completeness, 
these results can be translated into angular asymmetries, yielding A + = 0.53 ± 0.02 and 
A_ = -0.84 ±0.02. 

An alternative analysis, based on requiring two 6-tagged jets to further suppress the 
W+jets background and events with misidentified leptons, was used to measure the helicity 
fractions in the single-lepton channels as a cross-check. The observed cos 6* distributions 
were corrected by subtracting the expected background contributions and were unfolded 
using correction functions in an iterative method similar to the one used in the measurement 
of the angular asymmetries described in Section 4.2. The results were found to be in 
agreement with the single-lepton measurements presented in this paper. 

6.2 Constraints on the Wtb vertex structure 

Any deviation of Fq, Fl, (or A + and A J) from the Standard Model prediction could be 
caused by new physics contributing to the Wtb vertex. Such new interactions associated 
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Figure 6. Overview of the four measurements of the W boson helicity fractions and the combined 
values. The error bars correspond to the statistical and total uncertainties. 

with the top quark may exist at higher energies. New physics can be parameterized in 
terms of an effective Lagrangian [11] above the electroweak symmetry breaking scale of 
v = 246 GeV. After electroweak symmetry breaking, the Wtb Lagrangian [12, 55] is: 

C wtb = —jfr (V L P L + V R P R ) t W~ - -£=b (g L P L + g R P R ) t W~ + h.c. , (6.1) 



where Pl (Pr) is the left-handed (right-handed) chirality operator and 

v L = v tb + c£ 3+3) ^ , v R = 1 -cf;^, 9L = V2C 3 d ^, 9R = V2C 3 J V ^. (6.2) 

,(3,3+3) ^33* ^,33* , ^33 



The parameter A is the new physics scale and C\ ' , C^*, and C'^y are the 

effective operator coefficients [13, 55]. The anomalous couplings V R , <7l> 9r, generated by 
dimension-six operators, are absent in the Standard Model at tree level, while the coupling 
V t b receives a correction from the operator 0^ 3+3 ^ . 

Limits on anomalous couplings (V R , g\, and g R ) were obtained from the combined 
measurement of the W boson helicity fractions by exploiting their dependence on these 
couplings, as implemented in the TopFit program [10, 56] and normalizing to Vl = 1- 
The allowed regions of 9r) are shown in Figure 7, assuming V R = 0. The upper 
disconnected region in the plot shows a laige-g R second solution to the quadratic equation 
relating the observables to the anomalous couplings. However, this region is disfavored by 
the measured cross-section for single top production at the Tevatron [17, 57, 58]. 

In addition to this two-dimensional limit it is useful to set limits on single anomalous 
couplings, taking only one of them non-zero at a time. These are, at 95% confidence level, 
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Figure 7. Allowed regions at 68% and 95% confidence level (CL) for the Wtb anomalous couplings 
5l and <?r. In the Standard Model, the anomalous couplings vanish at tree level [59]. 



Re (C 33 ) 

Re(Vk) G [-0.20,0.23] -> G [-6.7,7.8] TeV~ 2 



A 2 



Re( 5L ) G [-0.14,0.11] -> G [-1.6,1.2] TeV^ 2 



Re(<7 33 i 

Re (<7r) G [-0.08,0.04] -> G [-1.0,0.5] TeV" 2 . 

The considered W boson helicity observables also allow a second region for gn when the 
remaining anomalous couplings vanish: Re (<7r) G [0.75,0.80] at 95% confidence level. 
It should be noticed, however, that such large coupling values would imply a single top 
production cross-section value disfavored by the Tevatron measurements [17, 57, 58]. Using 
a Bayesian approach [60], the measurement of the W boson helicity fractions with 
fixed at zero, was translated into a 95% probability interval on Re(C 33 y)/A 2 , as proposed 
in Rcf.[13]. This interval was found to be [—0.9,2.3] TeV~ 2 . 

It can be seen that the limits on Cjjy (mediating the production of right-handed b- 
quarks in the top decay) are of the same order of magnitude as the limits on C 33 y (involving 
left-handed quarks). This reflects a good sensitivity to the effective operator corresponding 
to Cjjy, even if its contribution is suppressed by 1/A 2 instead of 1/A [61]. 

These limits are more stringent than those obtained by the D0 Collaboration [58, 
62] 4 . Indirect, model-dependent limits on the anomalous couplings have been inferred from 
measurements of radiative £>-meson decays, measurements of l?i?-mixing and electroweak 



4 The limits from the D0 Collaboration were derived assuming a massless 6-quark. 
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precision data [63-67], but these limits include assumptions on the absence of additional 
new physics effects that are not needed in this analysis. 

7 Conclusions 

A measurement of the polarization of the W bosons in top quark decays was presented, 
based on 1.04 fb _1 of data collected with the ATLAS detector in 2011. The single-lepton 
and dilepton decay topologies of top quark pairs were considered in the analysis. 

The helicity fractions obtained from a combination of template fits to the reconstructed 
cos 9* distributions and angular asymmetries calculated from the unfolded cos 9* distribu- 
tions are F = 0.67 ± 0.07, F L = 0.32 ± 0.04 and F R = 0.01 ± 0.05. These results are 
in agreement with NNLO QCD predictions and are more precise than previous results 
obtained by the CDF and DO Collaborations [6-8]. 

Limits on the Wtb vertex anomalous couplings were obtained from the combined results 
on the W boson helicity fractions. These results are consistent with the (V — A) structure 
of the Wtb vertex and improve on the previously obtained limits [58] . 
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